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Abstract 

In the Standard Model (and MSSM), renormalization effects on neutrino mixing 
are generally very small and the attractive fixed points are at vanishing neutrino 
mixing. However for multi-higgs extensions of the Standard Model, renormalization 
effects on neutrino mixing can be large and nontrivial fixed points are possible. 
Here we examine a simple two-higgs model. For two flavors, maximal mixing is an 
attractive infrared fixed point. For three flavors, the neutrino mass matrix evolves 
towards large off-diagonal elements at low energies. The experimentally suggested 
bimaximal neutrino mixing pattern is one possible attractive infrared fixed point. 
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1 Introduction 



Recent experiments JT|, [| have revealed important features of the neutrino mass matrix. 
It is now rather well-established that the neutrino masses are tiny, and that at least some 
of the mixing angles are large, or even maximal. Considerable effort has been devoted to 
a theoretical understanding of these features. 

The small value of neutrino masses is a natural feature of theories with new physics 
entering at a high energy scale, such as Grand Unified Theories (GUT). In such theories, 
neutrino masses enter the low energy effective Lagrangian as a dimension five operator || 
involving two fermion fields and two higgs fields. Thus the neutrino masses are suppressed 
from the charged fermion masses by a factor of the ratio of the weak scale to the scale 
of new physics. However these dimensional arguments do not explain why the neutrino 
mixing angles should be large. 

Radiative corrections renormalize all terms in the Lagrangian, including the neutrino 
mass. These effects may provide a dynamical explanation for the value of the top quark 
and Higgs masses (for a recent review, see e.g. M). For neutrinos, these corrections were 
worked out some time ago || |6| and have been studied extensively of late for the Standard 
Model (SM) and the minimal supersymmetric standard model (MSSM) (see e.g. [[7], [| ||). 
In these models renormalization effects on neutrino mixing are proportional to the charged 
lepton masses and hence are relatively small. Also, maximal neutrino mixing is a saddle 
point in the SM and MSSM ||, with the infrared attractive fixed point corresponding to 
vanishing neutrino mixing. Thus renormalization effects in the SM and MSSM can not 
readily explain the observed large neutrino mixing. 

Here we shall study neutrino mixing in an extension of the Standard Model to two 
Higgs doublet fields, $i and $2- It is customary to assume that the two Higg's transform 
differently under a discrete or continuous (e.g. symmetry chosen so that each 



type of charged fermion couples to only one Higgs doublet. This insures that there 
are no flavor changing neutral Higgs couplings in the dimension four terms. However 
this symmetry does not necessarily mean that there is only one type of dimension five 
neutrino mass term. In general, there are four ways to combine two Higgs fields and two 
neutrino fields 0]. Of the four dimension five operators, one involves two $i fields, one 
involves two $2 fields, and the other two depend on both fields. The symmetry used to 
eliminate flavor changing neutral currents among the charged leptons may be chosen such 
that almost any desired combination of these operators are allowed. Here we choose to 
study the case where only the latter two operators are relevant: 

where II is the left-handed lepton doublet, m and n the generation indices and i,j,k,l 
are SU(2) indices. K mn is symmetric under interchange of the generation indices m and 
n, while £ mn is antisymmetric. When the Higgs' fields acquire vacuum expectation values 
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($°) = Vi/ v^2 the neutrino mass term is 

M v = —-KV1V2 = —-Kf 2 sin 2/5. (2) 

where v 2 = vf + v% = (246) 2 GeV and tan/3 = Vijv\. Because the neutrino mass term does 
not explicitly depend on £, this combination of operators is relatively simple to study. 

Renormalization mixes the two operators so both must be simultaneously evolved to 
study neutrino mixing. In the basis where the charged lepton mass matrix is diagonal, 
the evolution equations are M 



^ = [-3s 2 + 2A 3 + 2A 4 + % Q ,-l( 2/ 2 — 



3g 2 + 2A 3 + 2A 4 + S]K a p - -(y + y a )n a p 
+2(y}-y 2 a)Ui3}, 
^ = [-9gl + 2X 3 -2\ i + S}^ + l(y 2 +yl)U 

+l(y}-y>a P \- (3) 

Here t = j^lnfi, g 2 is the SU(2) gauge coupling constant, and y a = ^/2m a /vi = 
\f2m ,/v cos/5 is the Yukawa coupling for charged leptons of type a (e, /i or r). The 
quantity S is defined as 

S = Tr[3YX + 3Y}Y d + YX)- (4) 
where Y a is the 3x3 Yukawa coupling matrix for charged fermions of type a. The Aj's 
are parameters from the higgs potential 

C2H = -^($I$ 1 ) 2 -^($^ 2 ) 2 -A 3 ($ t 1 $l)($^ 2 ) 

-A 4 ($ t 1 $ 2 )(^$ 1) - [y ($I$ 2 ) 2 + h.c] (5) 

These equations assume $1 couples to the charged leptons, but do not depend on which 
Higgs field couples to which type of charged quark. 

In the limit that £ — > 0, the evolution equation for k decouples and becomes similar 
to the SM case. We concentrate here on studying the new effects on neutrino mixing 
introduced by the £ parameter. To simplify the analysis we assume that all parameters 
are real. 



2 Two neutrino flavors. 

The evolution of the scale of the neutrino mass matrix can be separated out from the 
evolution of the dimensionless parameters. For dimensionless physical parameters in the 
two-flavor approximation we use the mixing angle 9 

tari 2# = 2/V , (6) 



2 



a ratio of neutrino mass eigenvalues, m,, defined as 

m 2 + mi 



m>2 — mi 



and the ratio of the upper off-diagonal elements of the two dimension five terms 

V = ^- (8) 

Since the diagonal elements of £ vanish, this parameterization is complete. The particular 
form of z is chosen so that the fixed points are finite and simple 0. In this parameteri- 
zation, 

m 2 - mi / z - cos29 sin 29 \ 

2 V sin2fl z + cos 26 J ^ ' 

After a little algebra, one finds that the evolution equations for the dimensionless 
parameters are 

| = f^-w 

^ = y 2 (i C os2^ 2 -l) -2^ sin 2 20) 
1 

— = ?/ 2 sin2#(-z + cos2^) (10) 
at 4 



where 

C = -6g 2 2 - 4A 4 + 2(y 2 T + y 2 ) 

= (11) 

These equations and the dimension five operators possess various symmetries. For 
example, under the transformation 

e^\-e (12) 

accompanied by an interchange of the mass eigenvalues 

z — > — z (m 2 <-> mi) (13) 

the diagonal elements of the neutrino mass matrix, Eq. (Q), are invariant while the off- 
diagonal elements change sign. However the overall sign of the off-diagonal elements is not 
a physical observable since it may be absorbed in the unphysical phases with a redefinition 
of the neutrino wave function {v^ u T ) — ► (i/^, —v T \ This transformation changes the sign 
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of £ MT and k^ t simultaneously, so 77 is invariant. Thus the neutrino mass matrix and the 
RGE evolution equations respect this symmetry. 

The running of rj is independent of the neutrino physical parameters z and 9. Thus 
we shall use rj to characterize the possible motion of z and 9. The fixed points of z and 
9 are obtained by setting their derivatives in Eq. flKp equal to zero and solving. The 
stabilities of each fixed point are obtained by finding the eigenvalues of the Jacobian at the 
fixed point (see e.g. [fL^Q. The results of this analysis for decreasing t (i.e. approaching 



the infrared) are in the Table (for increasing t the stabilities are reversed). As the Table 
shows, the location of the fixed points are independent of 77, but their stability does depend 
on rj. Of the five fixed points, only three are physically distinct since the symmetry of 
Eqs. ( |12D and (|H|) maps two into two others. The fixed points at 9 = 0, z = +1 and 
9 = n/2,z = — 1 correspond to no mixing and a massless muon- neutrino. The fixed 
points at 9 = 0, z = —1 and 9 = n/2,z = +1 correspond to no mixing and a massless 
tau-neutrino. The fixed point at 9 = ir/4,z = corresponds to maximal mixing with 
equal magnitude but opposite sign mass eigenvalues. Note that the SM corresponds to 
rj — > for which large mixing is a saddle- node |J. 

To determine which of the possibilities in the Table is realized, we must determine the 
behavior of 77. Setting = in Eq. ( JlQ| ) yields the fixed points of 77, t)±, as 



rj* ± = -(b±VWT3) 

C_ -6gl - 4A 4 + 2(y 2 + yj) 

2y 2 2(tf-yJ) ' 1 ) 

Thus rj evolution has two fixed points — one positive and one negative. The stability of 
the fixed points is easily found by plotting — ^ versus rj, as in Fig. (1). The minus sign 
is included in the derivative to give us the stability for t decreasing . We see that the 
positive fixed point, r/1, is a repellor while the negative fixed point, r/* , is an attractor. 
This is true, regardless of the value of the parameters C or y 2 . If the initial (high energy) 
value satisfies rj < rA, then rj evolves directly to the fixed point. If the initial (high 
energy) value satisfies r\ > 77I then rj initially evolves toward positive infinity. However 
note that evolution through 77 = +00 to the attractive fixed point at 77* is possible when 
the denominator of 77, k^ T) evolves through 0. 

From the Table, the attractive fixed point value of 77 will cause maximal neutrino 
mixing to be an attractive fixed point when 

V- < -\ (15) 

This corresponds to 

— cos a 5 < 1.75 (16) 

V VSM J 

,2 _ „,2 / _2 a „„J „,2 



where y = yg M / cos j3 and y$ u is the Standard Model value for the squared charged 



lepton Yukawa coupling differences (t/| m = 1.0 x 10 4 for y 2 — y 2l s and y| M = 3.5 x 10 ' 
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for y 2 — y 2 ). Our lack of knowledge of the Higg's parameters tan/? and A 4 (especially 
its sign) makes it impossible to definitively evaluate this condition. However it appears 
that maximal neutrino mixing is an attractive fixed point for an extremely wide range of 
parameters. 

The neutrino mixing and mass evolution for the case of negative rj is graphically 
displayed in Fig. (2). There the fixed points and the direction field are plotted for 
decreasing t. Starting from some initial point specified by the high energy theory, the 
evolution of z and 9 follow a trajectory on this graph and the plotted unit arrows show 
the directions tangent to this trajectory. The vector field plot does not explicitly depend 
on t or y 2 but does depend on the value of rj. However for \rj\ >> 1 the explicit r] 
dependence also cancels out of the vector field plot. Here we use a fixed point value 
77* = —900 which corresponds to g\ = 0.4, A4 = 1 and y 2 = 3.5 x 10 -3 . Except at the 
extreme upper and lower edges of the plot, the evolution is towards the large mixing fixed 
point. 

In the Standard Model, the small size of y 2 tends to suppress the neutrino parameter 
evolution. To compensate for this small factor, previous analyses 0, ||, |7], |8| have usually 
focused on the possibility of a large neutrino mass degeneracy at high energies, i.e. \z\ > 
1/y 2 . Let us briefly consider the evolutionary behavior of large \z\. For large \z\, z changes 
much faster than 9 does. For \z\ » \r)\, the evolution equations, Eq. fllOD , are the same 
as those of the Standard Model and the approximate invariant is ]9[] 

sin2 29 , 

—^C (17) 

where ( is a constant, independent of t and y 2 . Qualitatively, the running is away from 
the Standard Model repulsive fixed point and toward the Standard Model attractive fixed 
point (see Table). For small initial (high energy) 9 values, running causes z > and 9 to 
decrease until z ~ f] when 77 effects become relevant and then drive the mixing towards 
maximal. For large initial 9 values (9 near 7r/2), running causes z > to initially increase 
and 9 to initially decrease, until 9 evolves through maximal mixing to 9 < 7r/4 where z 
starts to decrease. When z has decreased enough such that \z\ ~ \rj\ then rj effects become 
relevant and drive the mixing back toward maximal. Thus it is only when \rj\ > \z\ that 
the attractive fixed point nature of maximal mixing is apparent. Sufficient evolution will 
eventually produce \rj\ » \z\. 

For large \r)\, the z dependence of % in Eq. fllPf) disappears and evolution of the 
mixing angle simplifies considerably. Then the renormalization equations may be solved 
to give 

tan 29 min « tan 2# max Exp[2 I ™" dty 2 r/} (18) 

which describes the evolution of the mixing angle. Here t max and t m i n are the high energy 
and low energy values of t. This equation clearly shows that running does not change 
the sign of tan 29, but just increases its magnitude, i.e. 9 — > tt/4. From Eq. ([jj§), the z 
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evolution can be obtained using the large 77 approximation invariant 



—^Z^X (19) 
cos if 

where x is a constant, independent of 77, and so only slowly varying with t and y 2 . 

Using Eq. (pT8|) we can find the approximate conditions for significant evolution toward 
the maximal mixing fixed point. A change in tan2# greater than a factor of 10 requires 

,2/ 



6.1 < {y 2 (-v)) (20) 

where we have taken /i min = 10 2 GeV, /i max = 10 15 GeV and the () denotes a value averaged 
over the running between these points. This condition may be satisfied if \rj\ is large 
and/or tan f3 is large, since y 2 = y 2 SM / cos 2 (3. The validity of perturbation theory requires 
y 2 < 4n. It is interesting to note that if y 2 is small, i.e. |6#2 + 4A 4 | » y 2 and we take 
T] ~ r)*_, the attractive fixed point value (Eq. (|14])), then the dependence on the charged 
lepton Yukawa coupling disappears from Eq. (^0|) to give 



(A 4 ) > \ 



ln(10) fi/ 2 

: 6 \92 



2.4 (21) 



where we have used g 2 ~ 0.4. The validity of perturbation theory requires A4 < 8tt. 

The evolution of 77 may be solved for approximately. If we neglect the evolution of y 2 
and C, then in Eqs. ( |I~0D may be integrated to give 

(^-y(, -C) = EM2A< _ <){tmax _ tmm)] (22) 

Using this equation, we may estimate the condition for evolution to produce large values 
of rj. Let us take r] max w —1, r] min w ^- << —1. Then Eq. ( ^2|) reduces to 

|r/* I < Exp[2|/ 2 |77* |(t TOaa; - t min )}. (23) 

If ?/ 2 is small, i.e. |6g| + 4A 4 | >> y 2 , then Eq. (p3|) reduces to 

1.2 + 2A, 



In 



Z/ 2 



< [0.46 + 0.76A 4 ] (24) 



Unlike Eq. ([H]), this conditions still weakly depends on y 2 . The smaller y 2 is, the larger 
the magnitude of the attractive fixed point value 77* , and so more running is required 
which in turn requires a larger A4. 

In general, the amount the mixing angle evolves depends on two parameters that 
are not well known, A4 and tan/3. A complete analysis of all the constraints on these 
parameters, and how they relate to neutrino mixing angle evolution, is beyond the scope 
of the present paper. However we note that studies of the evolution of the two-Higgs 
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potential parameters (see e.g. |13| ) suggest that the Aj's are not generally driven to a 
particular fixed point, but that a range of values are consistent with experimental and 
theoretical constraints. Constraints on two-Higgs models also follow from studies of flavor 
changing neutral currents. In particular, the flavor changing neutral current decay b — > s 
+ 7 tightly constrains two-Higgs models where a different Higgs couples to the two types 
of quarks (type-II) (see e.g. [0]). However this constraint can be avoided by coupling $ 2 



to both quark types and $1 to the charged leptons. Then large tan (3 corresponds to large 
charged lepton Yukawa couplings, but approximately Standard Model values of the quark 
Yukawa couplings. Large values of y 2 would help produce large mixing angle evolution, 
but they also tend to enhance lepton universality violations. However a recent analysis 
of W and Z decay data [15j finds very weak limits on two-Higgs' models from lepton 
universality, with values of tan (3 up to and exceeding 100 allowed (this corresponds to y T 
exceeding 1). Thus large running of the neutrino mixing angle appears to be compatible 
with present constraints on two-Higgs models. 



3 Three neutrino flavors. 



The evolution of the physical mixing angles and mass ratios is quite complicated for three 
flavors. Here we do not give explicit formulae for their evolution but instead focus on the 
texture of the mass matrix produced by evolution. Then we plot renormalization effects 
obtained numerically for one particular scenario. 

The general nature of renormalization effects on the neutrino mass matrix in this model 
can be found by examining the relative running of diagonal to off-diagonal elements 




vl) 



(25) 



Here r\ a p is defined just as it was for two flavors, 



Vaf3 



(26) 



T] a /3 is antisymmetric so for three flavors there are now three independent 77's; rj efJi , f] er 
and ?7^ T . For two-flavor, this equation provides another way to understand the Table. 
For r] a /3 < — 1/4 for (3 > a, the coefficient of (Kpp/n a p) on the right-hand side is positive 
for (3 > a, and also for a > [3 since r\ a p is antisymmetric. Thus as t decreases the 
two diagonal elements decrease with respect to the off-diagonal elements. This leads to 
maximal mixing, in agreement with the Table. Conversely for 1/4 < 7] a p for (3 > a, 
both diagonal elements increase with respect to the off-diagonal elements, which drives 
the mixing away from maximal mixing and towards vanishing mixing. Standard Model 
like behaviour occurs for —1/4 < rj a p < 1/4. Then the cofficient of (itpp/Kap) on the 
right-hand side is positive for (3 > a but negative for a > (3. Thus as infrared energies 
are approached the lower (upper) diagonal element increases (decreases) with respect to 
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the off-diagonal term. Running can go through the point where the diagonal elements are 
equal, and maximal mixing occurs, and proceed on towards vanishing mixing. Thus the 
texture of the k matrix depends on the behaviour of the rj a ps. 

The dynamics of the r] a ps follows from Eqs. There we see that £ Q/ 3 and n a p only 
mix with each other so, aside from different Yukawa couplings, the evolution equation for 
each three-flavor r] a p is identical in form to the two-flavor evolution equation for 77 given 
in Eqs. fllOf) . Thus the general dynamics of the f] aj 3 for (3 > a are identical to what has 
been discussed previously and plotted in Fig. (1). They all evolve toward attractive fixed 
points, 77*^, given by the negative root of Eq. fll4]) . Thus evolution produces r? Q/ g's which, 
for (3 > a, are negative and readily less than —1/4 (see Table and Eqs. (|1^) and ([16])). 
Referring back to Eq. (|25|) , this implies that evolution toward low energies decreases the 
diagonal elements relative to the off-diagonal elements of the neutrino mass matrix. Then 
in the infrared limit the neutrino mass matrix approaches the form 



v a/3 



t*»fo% ierhtr 



(27) 



Here the ^*'s are the attractive fixed point values and £'s are their values at low energies. 

For comparison purposes, let's consider the relative running of two off-diagonal ele- 
ments 



d_ 

dt 



-(y 2 T - vl) + 2{yl - yl)^ - 2{y 2 T - y 2 e ) Ve 



:28i 



This equation shows that a cancellation occurs between the two terms proportional to rfs. 
Indeed, if we substitute in the fixed point value of the rfs for small Yukawa couplings, 
all large terms cancel out and the running is not significant. Thus the running of the 
off-diagonal terms is generally not a dominant process. 

The form of the neutrino mass matrix in Eq. (^) is the same as that produced in the 
popular Zee model |16|] . However note that the Zee texture is never exactly obtained from 
renormalization, i.e. the diagonal elements never completely vanish. How closely the Zee 
texture is approached depends on the amount of evolution, the size of the initial (high 
energy) parameters and the values of the two-higgs parameters tan/5 and A4. But this 
general texture is interesting because it can describe the experimentally desired bimaximal 



neutrino mixing ]T7| when n eil « n er >> k^ t . This particular relationship among the off- 



diagonal parameters is not caused by running, as Eq. (|28|) shows. However it may be 
realized if the £ aj a's at the high energy scale have the appropriate hierarchy of values. 

The possible size and nature of running are demonstrated in Figs. (3) and (4). In 
these calculations, 400 different neutrino mass matrices were generated by choosing the 
elements of K a p at the high energy scale to be real, random numbers evenly distributed 
between -1 and 1. The £, a ps at the high energy scale were chosen to be £ efl = 2600, 
£ er = 10 and £ MT = 0. The parameters £ efl and £ er were chosen to be larger than the 
K aj g's at high energies to insure that K efJ- and n eT would be the largest matrix elements at 
low energies. Additionally, these parameters were chosen because they satisfy £ e/J /£ e r ~ 
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Veti/Ver (actually their ratio is a little less than this to allow for running before the r] a ps 
reach their fixed points). The other relevant parameters were taken to be g\ = 0.4, A4 = 1, 
y T — 1 (which is equivalent to tan/3 = 100), and [— 3g| + 2A3 + 2A4 + S] = 7.8 (the last 
influences the overall neutrino mass scale evolution). These parameters were not evolved 
but instead held fixed because they mostly influence the overall amount of running so 
their precise values are not particularly important, and because the Aj's and tan/3 are for 
the most part unknown. 

Fig. (3) shows the mixing angles at the high and low energy scales. The mixing angles 
are defined following the convention of ref. [18 as 



4772 772 

and 



c3j 



to to 


— m\\ 


< 


co to 


— m 2 


\ m l 


— m\\ 


< 


CO to 





ATI 2 U 2 

*>'*>> -ffife (30) 

where U a i is the unitary mixing matrix. The neutrino mass eigenvalues are ordered such 
that 

(31) 
(32) 

so that solar neutrino oscillations are described by the m\ — m 2 mass squared differ- 
ence. Thus ip is the mixing angle relevant for atmospheric "z^-zV oscillations and u is 
the mixing angle relevant for solar neutrino oscillations. Fig. (3) clearly demonstrates 
that renormalization from high to low energies can simultaneously drive these angles to 
maximal. The running effects on U& are not shown, however for the same parameters 
the relevant mixing quantity, sin 2 20 = 4U 2 3 (1 — U 2 3 ), is driven from a range of values 
between and 1 to sin 2 20 < 0.05 at the low energy scale. This is consistent with the 
experimental limits from CHOOZ | I9| . Thus running can drive the neutrino mass matrix 



to a bimaximal form, consistent with present experiments. 

The approximate conditions for large mixing angle evolution were described in the 
two-flavor section in Eqs. (|20D , (pl|) and Eqs. (p3|) , fl24"D . In particular, these equations 
shows that large angle evolution results if A4 is large. Indeed, our numerical results 
remain essentially unchanged if the charged lepton Yukawa couplings are decreased to be 
their Standard Model values (tan/3 << 1) and A 4 is increased to A 4 > 16. This agrees 
with Eqs. ( pT|) and (0), except that the numerical results are somewhat better than 
the analytical estimates suggested by these equations because the numerical results used 
initial values of £ efl and £ er that were large in magnitude. Large mixing angle evolution 
could be maintained for smaller values of these initial £ eAt and £ eT by increasing the size 
of A4. 

Fig. (4) shows how running affects neutrino mass squared differences. The quantities 
m 2 — m 2 and m\ — m\ are each scaled by the largest neutrino mass squared and plotted 
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against each other. From a distribution of values at the high energy scale, evolution 
produces an inverted hierarchy of neutrino masses with, roughly, m\ « m\,m\. This is 



consistent with the results found in studies of bimaximal mixing in the Zee model [17 



4 Conclusions 

Neutrino masses enter into low energy physics as a dimension five operator, regardless of 
the details of the high energy theory. In multi-higgs extensions of the Standard Model, 
there are generally more than one dimension five lepton-higgs operator. Here we have 
worked in the two-higgs model, picked two of the dimension five operators (one a sym- 
metric mass term, the other an antisymmetric operator), and examined how mixing angle 
evolution is changed from the Standard Model. Renormalization mixes these different 
operators, thus neutrino mixing depends on how the two operators evolve. Because the 
operators' evolution have different dependences on a higgs potential parameter, the mix- 
ing angle evolution depends on this parameter also. This parameter can be much larger 
than the tiny Standard Model charged lepton Yukawa couplings, so the mixing angle evo- 
lution can be quite large. The size of angle evolution is also enhanced because multi-higgs 
models have larger charged lepton Yukawa couplings. 

For the model studied here, for a wide range of parameters, running tends to produce 
a neutrino mass matrix whose diagonal elements are small compared to the off-diagonal 
elements. For two neutrino flavors, this corresponds to maximal mixing as the attractive 
infrared fixed point. For three neutrino flavors, the neutrino mass matrix resembles that 



produced by the Zee model |16|]. It is known that this model can describe the experi- 
mentally suggested bimaximal mixing. Here we showed explicitly that bimaximal mixing 
can be an attractive infrared fixed point if the elements of the asymmetric dimension five 
operator are large and heirarchical. 

This model predicts that the neutrino mass spectrum is inverted, with m| << m|, m\. 
This possibility may be tested using matter effects , either in long baseline experiments 
or in atmospheric neutrino observations. The scenario described here may be distinguished 
from the Zee model by observing experimental processes such as double beta decay and 
CP violation. These quantities are highly suppressed in the standard Zee model, but are 
expected to only be moderately suppressed by renormalization effects. 
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Table. The stability of the fixed points in our two-Higgs model for t decreasing, i.e. 
the infrared limit. The notation is A = attractor, S = saddle, R = repellor. 
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Figure 1: Plot of — % versus rj, the ratio of off-diagonal terms between the two dimension 
five neutrino-higgs operators. We have taken g\ = 0.4, y 2 = 3.5 x 1CT 3 and A4 = 1. The 
solid circle and open circle denote the attractive and repulsive infrared fixed points. 
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Figure 2: Direction field for running of the physical neutrino parameters. Two neutrino 
flavors are assumed with z being the ratio of mass differences, 9 the mixing angle and 
rj = —900. The different infrared fixed points are shown with solid circles, open circles 
and grey square denoting attractors, repellors and saddle point. 
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Figure 3: The plots show how evolution affects the solar sin 2 2cj and atmospheric sin 2 2ip 
neutrino mixing parameters. Random neutrino mass matrices at the high energy scale 
(fi = 10 15 Gev) are evolved down to the low energy scale (/i = 10 2 Gev). The fixed 
parameters are g\ = 0.4, A4 = 1, y T = 1, £ e/i = 2600, £ eT - = 10 and ^ T = 0. 
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Figure 4: The plots show how evolution affects the neutrino mass squared differences. 
Random neutrino mass matrices at the high energy scale (/z = 10 15 Gev) are evolved down 
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the previous figure. The mass squared differences are scaled by the largest neutrino mass 
squared. 
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